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A weighted average weak mixing angle Ow derived from relatively low experiments is compared 
with the Standard Model prediction obtained from precision measurements. The approximate 1.8 
sigma discrepancy is fit with an intermediate mass (~ 10 — 35 GeV) “dark” Z boson Zd, corre¬ 
sponding to a U{l)d gauge symmetry of hidden dark matter, which couples to our world via kinetic 
and Z-Zd mass mixing. Constraints on such a scenario are obtained from precision electroweak 
bounds and searches for the rare Higgs decays H —>■ ZZd —>■ 4 charged leptons at the LHC. The 
sensitivity of future anticipated low measurements of sva? 6w{Q^) to intermediate mass Zd is 
also illustrated. This dark Z scenario can provide interesting concomitant signals in low energy 
parity violating measurements and rare Higgs decays at the LHC, over the next few years. 


Discovery of what appears to be a fundamental Higgs 
scalar mm completes the basic Standard Model (SM) 
particle spectrum. In addition, comparing precision fine 
structure constant a, Fermi constant Gp, and Z boson 
mass (mz) values at the quantum loop level, employing 
the Higgs mass mn = 125 GeV and top quark mass 
rtit = 173.3(8) GeV gives the indirect SM weak mixing 
angle prediction mm 

sin^ = 0.23124(12) SM prediction, (1) 

where the modified minimal subtraction (MS) definition 
at scale fj- = 'm-z for the renormalized weak mixing angle 
9w has been employed [5]. The existing error in Eq. Q 
stems from m^, higher order loops (that overall double 
the error), and hadronic uncertainties, all of which are 
expected to be further reduced. That prediction agrees 
remarkably well with the average value j3] of the more 
direct Z pole measurements 0111 

sin^ 6 *u/(m. 2 )^ = 0.23125(16) Z pole average. (2) 

A comparison of these distinct precision methods severely 
constrains “new physics” extensions of the SM [3]. 

In contrast, low determinations of the weak mixing 
angle (for a review, see Ref. [3]) currently allow consider¬ 
able room for certain types of new physics, particularly 
Z' bosons (for earlier work along these lines, see for ex¬ 
ample Refs. [ 8 HTT]L Indeed, the 3 most precise measure¬ 
ments at lower rrv^ extrapolated, for comparison, 

to an MS scale /i = rriz give a somewhat disparate range 
of values [3] 

sin2 6 »w(mz)j^ = 0.2283(20) APV, (3) 
sin^ 0w{mz)^ = 0.2329(13) Moller E158, (4) 

sin2 0vv(TOz)Ms = 0-2356(16) NuTeV (5) 

from the measurements in Gs atomic parity violation 
(APV) at (Q) = 2.4 MeV [T3HT3] . SLAG Moller scatter¬ 
ing experiment E158 at (Q) = 160 MeV [16], and Permi- 
lab neutrino deep inelastic scattering (DIS) experiment 
NuTeV at (Q) « 5 GeV [17]. 
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FIG. 1. Current measurements of the weak mixing angle at 
various Q (6] [T] 1131419) and future prospects |20H24j . The 
black curve represents the expected SM prediction for the run¬ 
ning of sin^ 9w with Q [5] . Current measurements are given 
as black points with existing error bars. The red “Antici¬ 
pated sensitivities” are meant only to illustrate the possible 
uncertainties potentially obtainable from experiments under 
analysis and proposed. 


These measurements are illustrated in Fig. [^ af¬ 
ter evolving back to their experimental Q values. 
There, we also show other less precise determinations of 
sin^ dwiQ^) (JLAB Qweak first result [TS] and JLAB 
PVDIS [H]) as well as the very accurate Z pole val¬ 
ues mm, future sensitivities (Ra+ APV [20l |2T] , JLAB 
Moller [^, MESA P2 [^, JLAB DIS experiment 
SOLID [23]), and the predicted SM running curve for 
comparison. Note that the Qweak result in our figures 
corresponds to only about 4% of their total collected 
data. Their statistical uncertainty may be significantly 
reduced in the near future making them the expectedly 
best low determination. We return to this point later. 
Note, also, that the factor of 5 improvement envisioned 
for APV using single ionized Ra+ trapped atoms as orig¬ 
inally suggested in Ref. [33] , although extremely well mo¬ 
tivated, is still in a development stage [26]. The potential 
polarized electron scattering asymmetry improvements 
are currently on a more definite footing. 
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The weighted average from Eqs. (§-§ 
sin^ = 0.2328(9) low average (6) 


is roughly 1.8 sigma higher than the SM prediction in 

Eq. 0 


Asin^ Ow - 0.0016(9) 


(7) 


and gives about the same deviation relative to Eq. 0. 

Of course, there are still outstanding issues regarding 
atomic parity violation theory [27H29j that warrant fur¬ 
ther scrutiny. In addition, NuTeV hadronic effects [30] 
and radiative corrections [311 Ea could shift the average 
somewhat 3 . However, here, we take the current aver¬ 
age in Eq. (61 at face value and examine its consequences 
for an intermediate mass dark Z [Z^) with mz^ 

10 — 35 GeV (the intermediate mass range bounded from 
below by the onset of severe constraints from low energy 
measurements and from above by rriH—mz) and coupling 
to the SM particles via kinetic and Z-Z^ mass matrix 
mixing. Although the current 1.8 sigma discrepancy is far 
from compelling evidence for “new physics”, it does merit 
watching as low measurements of siv? 9 y/{Q^) along 
with independent constraints on Zd mixing improve. 

We start our discussion of intermediate mass Zd by 
briefly recalling its basic features. That scenario assumes 
a U(l)d gauge symmetry associated with a hidden dark 
sector. Its gauge boson, Zd, couples to our world (SM) 
via kinetic mixing, parametrized by £, and Z-Zd mass 
matrix mixing, parametrized by ez = {inzj^/'mz)5 [33] 
Actually, for an intermediate mass Zd, the combination 


5' ~ 5 -I- £ tan Ow 

mz 


( 8 ) 


proves important, as it governs the induced weak neutral 
current interactions of Zd (throughout our discussion, we 
ignore higher order corrections in e and <5). It means the 
5 is replaced by the more general 6' of Eq. 0 for an 
intermediate mass Zd- For the usually considered case of 
mzj mz, the second term in Eq. 0 [33] is generally 
negligible and 5' ^ 6 becomes a good approximation, 
but here it is retained. Depending on the relative sign of 
5 and e, the Z-Zd mass mixing or 5' might increase or 
decrease as mz^ increases. 

As a result of mixing, Zd couples to the SM via [33] 


Aint = (-e£Jr - 


5' 

2cos9w mz ^ 


+ (9) 


where the ellipsis represents other induced Zd interac¬ 
tions such as the HZZd coupling [331 EH ES] that we 
subsequently employ. As a consequence of Eq. 0, weak 
neutral current SM amplitudes at low momentum 


^ We note that a new Higgs doublet charged under U{l)d, assumed 
in typical models of Z-Zd mass mixing discussed in Ref. [33], can 
also lead to non-zero kinetic mixing, via loop effects. 



FIG. 2. BR(7y —>■ ZZd)l5''^ with For the most part 

{mzd ^ 30 GeV), the branching ratio into ZZd is almost 
independent of mz^- BR{H —>■ ZZd) ~ (16 — 18) S'^. 


transfer are rescaled by pd (that is pdGp instead oi Gp) 
and the SM weak mixing angle sin^ 9w{Q^)sm is replaced 
by Kd sin^ 9w{Q^)sm [^EZIEE] with 


Pd — 1 + 5 '' 




■ m 


Zd 


and 


( 10 ) 


Kd = l- eb -cot 9w yyt- — 

mzd Q 


( 11 ) 


''Zd 


The above yields a low <C shift 


A sin^ 9w — —£5' cos 9w sin 9w 


mZd 
~ -0.42 £(5' 


mz 

Tnzd 


( 12 ) 


Note that the effect of pd in Eq. (10) on sin^ 9w{Q^) is 
process dependent. Its largest effect is on the NuTeV re¬ 
sult of Eq. 0 , where an upward shift in the experimental 
sin^ 9w{m^-^ of 5'"^ is induced if Ri, (the ratio of neu¬ 
tral current to charged current neutrino cross sections) is 
employed [33 Ea, and (5'^/2 if the Paschos-Wolfenstein 
relation [33] is used. Overall, pd has little effect on the 
weighted average in Eq. 0. Nevertheless, including the 
effect of Pd in future more prec ise studies is warranted. 

As can be seen from Eq. (12), the value of sin^ OwiQ"^) 
in our framework depends on mz^, £, and 6'. Let us 
then consider next the current constraints on the latter 
two quantities over the mZd range of interest here. 

Recently, the ATLAS collaboration at the LHC has 
reported results for the rare Higgs decay H —>• ZZd 

with £i ^2 = e,/i [40]. Assuming Z-Zd mass 
mixing parametrized by S' and a dominantly SM-like 
Higgs boson of 125 GeV, one can show [33] that this 
decay has a branching ratio (roughly including Zd phase 
space effects [36]) 


BR(iJ ^ ZZd) « (16 - 18) S'' 


(13) 
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FIG. 3. Effective weak mixing angle running as a function of sh ift ( the blue band) due to an intermediate mass Zd for (a) 
mzd = 15 GeV and (b) =25 GeV for 1 sigma fit to eS' in Eq. (121. The lightly shaded area in each band corresponds to 


choice of parameters that is in some tension with precision constraints (see text for more details). 


which is further reduced by Z and leptonic branching 
ratios. The on-shell branching ratio is given by [33113S] 


BR{H -)■ ZZd) = 




gmz 



with A(a;, y, z) = x'^ + y'^ + z'^ — 2xy — 2yz — 2zx and 
rjj(125 GeV) ~ 4.1 MeV jH], which shows a rather mz^ 
independent value over most of the mass range (Fig. |^, 
resulting in Eq. (13). 


The ATLAS bounds translate into constraints on S' 
as a function of mz^, but depend on the branching ra¬ 
tio for Zd i+r. For BR{Zd 2i) = RR{Zd 
2e) +BR{Zd 2y) « 0.3 [13], one finds (at 2 sigma) the 
nearly constant bound |5'| < 0.02, over the range of mz^ 
considered in our work. Here we note that in the pres¬ 
ence of allowed dark decay channels (that is, decay into 
invisible particles), BR{Zd —>■ 2() can be much smaller 
than 0.3, which would weaken the constraint on S'. 

The best current bounds on e for the relevant mass 
range are given by the precision electroweak constraints, 
along with the non-continuous bounds from the e“'"e“ —>■ 
hadron cross-section measurements at various experi¬ 
ments [43]. The Drell-Yan dilepton resonance searches 
at the LHC experiments (such as in Refs. [HI US]) have 
the potential to give a better bound than precision elec¬ 
troweak constraints |4f)j . When combined with bounds 
on e from precision measurements and production con¬ 
straints 031117], one finds |e| < 0.03, for kinetic mixing 
alone. However, in our scenario, where a separate source 
of mass mixing is also considered [33] , that bound can be 
somewhat relaxed, via partial cancellation with S' depen¬ 
dent contributions to the Z-Zd mixing angle |33j , roughly 
yielding |e| < 0.04. (See also Refs. 07] US] for less severe 


bounds on e from a recasting of a CMS analysis of Run 
1 data, sensitive to H —>■ ZZd-) 

Given the above discussion, a simple combination of 
the upper bounds on e and S' suggests 


|e(5'| < 0.0008. 


(15) 


We use the above bound as a rough guide for the allowed 
region of parameter space in our discussion below. 

For a given mz^^ a neg^ive eS' in Eq. ([T^ will shift 


the SM prediction in Eq. ([^ towards the low experi¬ 
mental sin^ 0w(m2)j4g weighted average in Eq. (j^. That 
effect is illustrated in Fig.|^(a), where for mz^ = 15 GeV 
the blue band corresponds to a 1-cr fit to Eq. 0 or 
—0.0010 < eS' < —0.0003. A similar l-cr band is pre¬ 
sented in Fig. 1^ (b) for mz^ = 25 GeV with —0.0016 < 
eS' < —0.0005. In each case, the lighter shaded upper 
part of the band corresponds to |e(5'| > 0.0008 which 
is in some tension with constraints from precision mea¬ 
surements and the rare Higgs decay search by ATLAS, as 
explained above. Future improved sensitivity at the LHC 
should cover most of the bands in Figs. (a) and (b). For 
other rnz^ values, the 1-cr bands are about the same as 
our Fig. [^ representative examples; however, for larger 
mzd >25 GeV, the darker parts of the bands allowed 
by current constraints narrow. This can be seen from a 
comparison of Figs. 3 (a) and (b) that shows how smaller 
values of mzd can accommodate a shift in sir? 9w{Q'^) 
more easily, over the currently allowed parameter space 
[as suggested by the mz^ dependence in Eq.(12)]. 

In the case of low determinations of sir? 
the Qweak polarized ep asymmetry experiment at JLAB, 
which measures weak nuclear charge of proton (Q^gak)’ 
is expected to reach an uncertainty of ±0.0007 after all 
existing data are analyzed in the near future. This would 
reduce the uncertainty on the weighted average in Eq. 0 
to ±0.00055 and, assuming the same central value as the 
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current published result, could yield a ^ 3 cr deviation 
from the SM result in Eq. Q. It will be interesting 
to watch that outcome. We note that the weak mix¬ 
ing angle extracted from the Qweak experiment will ex¬ 
hibit some dependence on nucleon form factors including 
strangeness matrix element effects [49l [50]. For that rea¬ 
son, lattice gauge theory improvements in those hadronic 
matrix elements are strongly warranted. 

Future experiments, primarily polarized e e Moller 
scattering at JLAB and polarized ep scattering (P2) at 
MESA in Mainz, are expected collectively to further re¬ 
duce the weighted average uncertainty on sin^ 
at low below ±0.0002, becoming competitive with Z 
pole measurements. Together, low precision studies 
combined with improved H ^ ZZd searches at the LHC 
will squeeze the intermediate mass Zd scenario with some 
possibility of uncovering its existence. 

The intermediate mass Zd is an interesting viable al¬ 
ternative to the “light” dark photon often considered in 
the literature m- In addition to the parity violation at 
low that we have explored, it can give rise to poten¬ 
tial signals at the LHC, both in direct Drell-Yan produc¬ 


tion pp ZdX or as a, final state in rare Higgs decays. 
Besides the H —)■ Z Zd mode that we have discussed, 
searching for the mode H —>• ZdZd, mediated by Higgs- 
dark Higgs mixing [33], is well motivated. In fact, we 
note that the ATLAS 8 TeV search for H ^ ZdZd has 
two interesting but tentative candidate events (each at 
1.7 cr), roughly in the mass range ^ 20 — 25 GeV [30] . 
Further data from Run 2 at the LHC will be needed to 
clarify whether these events could be identified as inter¬ 
mediate mass Zd states that connect our world to an as 
yet unknown dark sector of Nature. Such a discovery 
would certainly revolutionize elementary particle physics 
and perhaps provide a new window into the world of dark 
matter. 
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